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We show how to ehminate the first-order efltects of the spin-orbit interaction in the performance 
of a two-qubit quantum gate. Our procedure involves tailoring the time dependence of the cou- 
pling between neighboring spins. We derive an effective Hamiltonian which permits a systematic 
analysis of this tailoring. Time-symmetric pulsing of the coupling automatically eliminates several 
undesirable terms in this Hamiltonian. Well chosen pulse shapes can produce an effectively isotropic 
exchange gate, which can be used in universal quantum computation with appropriate coding. 



PACS: 03.67.Lx, 71.70.Ej, 85.35.Be 



The exchange interaction between spins is a promis- 
ing physical resource for constructing two-qubit quan- 
tum gates in quantum computers In the ide- 
ahzed case of vanishing spin-orbit coupling, this in- 
teraction is isotropic, and any Hamiltonian describing 
time-dependent exchange between two spin-1/2 qubits, 
ffo(0 — J{t)^i • §2, commutes with itself at different 
times. Thus, the resulting quantum gate depends on J(t) 
only through its time integral — a convenient simplifica- 
tion, particularly because, when carrying out quantum 
gates, the exchange interaction should be pulsed adia- 
batically on time scales longer than fi/AE, where IS.E 
is a typical level spacing associated with the internal de- 
grees of freedom of the qubits |^. In addition, isotropic 
exchange alone has been shown to be sufficient for uni- 
versal quantum computation, provided the logical qubits 
of the computer are properly encoded [HQ]. 

Given the potential advantages of isotropic exchange 
for quantum gates, it is important to understand the ef- 
fect of the inevitable anisotropic corrections due to spin- 
orbit coupling. When these corrections are included, the 
Hamiltonian describing time-dependent exchange is 



H{t) = J{t){Si-?,2+A{t)). 



where 



A{t) ^ I3{t) ■ {Si X S2) + Si • ir(t) ■ S2 



(1) 



(2) 



Here (3{t) is the Dzyaloshinski-Moriya vector, which is 
first order in spin-orbit coupling, and W{t) is a symmet- 
ric tensor which is second order in spin-orbit coupling 
1^. Although these corrections may be small, they will, 
in general, not be zero unless forbidden by symmetry. For 
example, Kavokin has recently estimated that (3{t) can 
be as large as 0.01 for coupled quantum dots in GaAs 

In this Letter we construct the quantum gates pro- 
duced by pulsing H{t). This is nontrivial because H{t) 
typically does not commute with itself at different times. 
We represent the resulting gates using an effective Hamil- 
tonian H{t), which we derive perturbatively in powers of 



the spin-orbit coupling. H{t) is simple to work with be- 
cause it does commute with itself at different times. As 
an application of this effective Hamiltonian, we use it 
to tailor pulse forms that effectively eliminate any first- 
order anisotropic corrections. 

The quantum gate obtained by pulsing a particular 
H[t) is found by solving the time-dependent Schrodinger 
equation i-^\^i{t)) ^ H{t)\'ii{t)) where \^{t)) is the state 
vector describing the two spin-1/2 qubits (here, and in 
what follows, h ^ 1). In general this problem cannot be 
solved analytically. However, since we expect spin-orbit 
coupling to be small, it is natural to attempt a perturba- 
tive solution in powers of /3(t) and Tr{t) . To do this it is 
first necessary to solve the unperturbed {A{t) = 0) prob- 
lem exactly. This corresponds to pulsing the isotropic 
exchange interaction, for which the unitary time evolu- 
tion operator at time t is 



Uoit) = Texp 



J{t')Si ■ S2 dt' 



exp {—ix{t)Si ■ S2) 



where 



x{t) 



J(t')dt'. 



(3) 



(4) 



Here T is the usual Dyson time ordering, and the second 
equality in (^) follows from the fact that isotropic ex- 
change commutes with itself at different times. The un- 
perturbed quantum gate produced by a full pulse is then 
Uo{t — *■ 00) = exp(—iASi • S2) where A = J[t)dt is 
the pulse strength. This is a well studied class of quan- 
tum gates 0]. For A = tt the result is a simple swap, 
and for A = Tr/2 it is a 'square root of swap' which, in 
conjunction with single-qubit rotations, can be used to 
construct a controlled- not gate |^]. 

We now consider the effect of the anisotropic correc- 
tions A{t). Given the evolution operator for the un- 
perturbed system, it is possible to recast the problem 



1 



in the interaction picture by introducing the state vec- 
tor \^i{t)) = UQ(t)\^{t)) which satisfies the Schrodinger 
equation = J(i)^/(i)|«'/(i)) where Ai{t) = 

U^{t)Ait)Uo{t). A formal expression for the unitary op- 
erator describing a full pulse in this picture is then 



Ui = Texp —i 



J{t)Ai{t)dt 



(5) 



Expanding the exponential in generates the standard 
time-dependent perturbation theory expansion for [// in 
powers of Ai{t). Returning to the Schrodinger picture, 
the unitary operator describing the full quantum gate is 
U = exp(-aSi • S2)C//. 

Rather than simply carrying out the perturbation ex- 
pansion for Uj, it is useful to parametrize the resulting 
quantum gate in terms of an effective Hamiltonian of the 
form 



H{t) = J(i)(Si -82 + ^), 



where the time dependence of J{t) is the same as in 
and A is independent of time. Unlike H{t), the effec- 
tive Hamiltonian H{t) commutes with itself at different 
times. Thus, after a full pulse, H{t) yields the quantum 
gate U = exp (— zA(Si • S2 + A)). Our goal is then to 
find the operator A for which U is equal to the quantum 
gate produced by a full pulse of H{t). 

Because H{t) is traceless at all times t, the correspond- 
ing unitary time evolution operator has determinant 1, 
i.e., U e SU[4]. Requiring that our effective Hamilto- 
nian produce the same quantum gate then implies that 
A must also be a traceless Hermitian operator. The most 
general such operator acting on the Hilbert space of two 
qubits can be written 



^ = /3-(Si xS2) + Si-r-S2 



I . (Si - S2) + 1 • (Si 



where r is a symmetric tensor. This can be seen by 
noting that A is indeed traceless and Hermitian, and has 
15 independent real valued parameters, the number of 
degrees of freedom for a 4x4 traceless Hermitian matrix. 

Before proceeding it is instructive to classify the terms 
in A according to their symmetry properties under in- 
version (Si S2) and time reversal (S^ — > — Sj). Under 
inversion (3 and a change sign, while IT and Ji do not. 
Since /3(t) also changes sign under inversion this implies 
that (3 and a. are first order in spin-orbit coupling, while 
IT and Jl are second order. Under time reversal a and 
71 change sign, while /3 and IT are unaffected. We there- 
fore expect that for time-reversal symmetric pulses, i.e., 
pulses for which H{t[) — t) = H{t) (where is the center 
of the pulse), a and Jl will vanish. 

To determine A for a given pulse we note that the re- 
quirement that U = U implies 



(6) Texp 



J{t)Ai{t)dt 

= Texp(^-iJ J{t)Ai{t)dtj , (8) 

where Ai{t) = UQ{t)AUo{t). Expanding both sides of 
(^) to a given order in spin-orbit coupling and equating 
matrix elements yields a set of 15 independent equations. 
These equations can then be solved for the parameters 
in A in terms of J{t), f3{t) and W{t). 

We have carried out this calculation to obtain the fol- 
lowing expressions valid to second order in spin-orbit cou- 
pling (i.e., second order in /3{t) and first order in W{t)), 
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2sin(A/2) 



(3{t)sm[x{t)-^]j{t)dt, (9) 



S2), (7) /3 = 



J-^jy^ cos {.(t)-^)jm, (10) 



M = ^ / Jiti)dh I J{t2)dt2 {{I3{h) X /3(t2)) cos(a;(ti) - x{t2)) + 2(a x /3) sin{x[h) - x(t2))) , 



(11) 



and 



where 



Tab=jJ Tab{t)J{t)dt + ^ / J{h)dh I J{t2)dt2lab{hM) Sm {x{ti) - x{t2)) 



labihM) = 2 (f3{h) ■ /3(i2) - /3' - a^) 5ab - {Paitl)Pbit2) + Pait2)Pbitl) " 2/3^/3, - 2aaab) 



(12) 



(13) 



The criterion for the validity of these expressions is erator U. It is then apparent that, for any finite (3{t) 
that I A/3 1, |Aa| <C 1, where the factor of A is included and ir(t), our expansion breaks down when A 27rn for 
because it is the product A^ that enters the unitary op- n — ±1, ±2, • • •, because sin(A/2) ^ at these points. 
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However, for A — > 0, while a. and (3 may diverge, Acf and 
A/3 will always remain finite, and so, provided (3{t) and 
ir(t) are small, our expansion remains valid in this limit 

As expected from symmetry considerations, we find 
that /3 and a. are first order in spin-orbit coupling, while 
IT and 'p, are second order. It is also readily verified that 
for a time-reversal symmetric pulse the integrals and 
( pi] ) for a and /I vanish. Thus these non time-reversal 
symmetric terms are only generated by pulses that are 
themselves not time-reversal symmetric. 

Given the possibility of using the exchange interac- 
tion alone to perform universal quantum computation 
1^,0 , which depends crucially on the interaction being as 
close to isotropic as possible, a natural questions arises: 
Is it possible to ameliorate the effect of spin-orbit in- 
duced anisotropy on exchange-based quantum gates? We 
show below that the answer is yes — by carefully shap- 
ing pulses, it is possible to effectively eliminate the first- 
order anisotropy terms leaving only a residual second- 
order anisotropy. 

There are two first-order terms in H{t)^ a and (3. We 
have already seen how to eliminate a. By choosing a 
time-reversal symmetric pulse both a and Jl will vanish 
from H{t). Although (3 cannot similarly be eliminated, 
for appropriate pulse forms it can be effectively elimi- 
nated by performing a local rotation in spin space. 

Let S2 = IR • S2 where IR is a rotation matrix con- 
structed to eliminate (3 from H{t) so that 



H{t)^j{t) fsi -sJj + Si -ir'-sO , 



(14) 



where IT is a symmetric tensor. The precise form of this 
rotation depends on both (3 and IT and cannot be ex- 
pressed simply. However, up to second order in /3, it is 
given by 



IR-afc 



{p'5ab-PaPb)l'^ + O{0'), (15) 



and this is sufficient for our purpose of eliminating first- 
order anisotropy. Using (^Sj) one finds the residual 
anisotropy in (n4) is, up to second order in /3, 



r 



ah 



(/5'<5ab - /3,/5,)/2 + 0(/). 



(16) 



Thus, in this rotated coordinate system the first-order 
anisotropy vanishes and all corrections to the isotropic 
exchange interaction are second order in spin-orbit cou- 
pling. _ _ 

The ability to eliminate (3 from H{t) by simply ro- 
tating one qubit with respect to the other indicates a 
procedure for eliminating the first-order effects of spin- 
orbit coupling in any quantum computer that uses tun- 
able exchange for quantum gates. Suppose that symmet- 
ric pulses are used, so that a = 0, and pulse forms are 



chosen so that (3 is the same for all pulse strengths A. 
Then, if the qubits in the computer form a linear array, 
or any arrangement for which there are no closed loops of 
qubits connected by two-qubit gates, it will be possible to 
define a local spin-space coordinate system in which the 
effective interaction between any two neighboring qubits 
has the form (p^. While this procedure does not com- 
pletely eliminate the anisotropy, it does reduce it from 
an effect that is first order in spin-orbit coupling to one 
that is second order. 

To demonstrate how (|l^) can be used to tailor pulse 
shapes that lead to the same (3 for all pulse strengths A, 
consider the family of pulses 



j(i;A) = Jo(A) sech2(2t/T(A)), 



(17) 



where Jq{X) and r(A) are, respectively, the pulse height 
and width, and the pulse strength is A = J[t\ A) = 
Jo(A)r(A). To evaluate ( p^ ) it is also necessary to know 
the time dependence of /3(t). Determining the precise 
form of this dependence will require a detailed micro- 
scopic study of the specific realization of the exchange 
interaction being considered. Here we take, as the sim- 
plest possible illustrative model, a linear dependence on 
J(i;A), 



/3(<)=/3iJ(t;A), 



(18) 



for which the integral (^0|) can be performed analytically, 
with the result 



(3 = f3,^^{2-\cot{\/2)) 



(19) 



Also, because these pulses are time-reversal symmetric, 
(^) gives OL — 0. 

Equation (|9|) can be used to exploit the freedom to 
choose Jo (A) and r(A), while keeping Jo(A)r(A) — A, 
to shape pulses that keep (3 fixed for different pulse 
strengths. For example, if the pulse parameters for X — tt 
(swap) are fixed to be Jo(7r) and T(7r), then, for general 
A, one should take 



Jo(A) = Jo(^)- 



2A2 



and 



2- Acot(A/2)■ 



^(A)=T(7^)-(2-Acot(A/2)). 



(20) 



(21) 



These pulse forms are shown in Fig. ^ for various val- 
ues of pulse strength A. Note that as A increases, the 
pulse height decreases. This is because ^ becomes in- 
creasingly sensitive to f3{t) with increasing A until, in 
the limit A 27r, the pulse height must go to zero if (3 
is to be kept constant. Although our perturbation ex- 
pansion for A breaks down as A — > 27r, for this example 
the pulse heights are chosen so that the parameters in A 
remain small, and we are always within the perturbative 
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regime. The pulse forms defined by ( pO| ) and (21) are 
therefore valid, even in this singular limit. Of course, in 
practice, pulses near X — 2n will be problematic because 
of the diverging pulse length. 



1.25 



to eliminate (3 we find, using (|16|), that the residual 
anisotropy in H(t) is 
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FIG. 1. Pulse forms tailored to produce the same (3 for 
different pulse strengths A for the example described in the 
text. Of the pulses shown, the narrowest with the highest 
peak is for A = 7r/4. A then increases in increments of 7r/4 as 
the peak height decreases until, for the widest pulse with the 
lowest peak, A = 77r/4. As A ^ 27r the pulse height goes to 
zero. 

Once the first-order corrections to H{t) are eliminated, 
the residual second-order anisotropy can be found by first 
evaluatin g (p^ ) and then performing the local rotation to 
eliminate /3. As a specific example, consider the special 
case for which the form of the pulsed Hamiltonian is 

H{t) = J{t)Si-mit)-S2, (22) 

where IR(t) is a time-dependent rotation matrix. Such 
rotated exchange is, in fact, precisely the form of 
anisotropy found microscopically when spin-orbit correc- 
tions are included in the usual Hubbard model treatment 
of superexchange |liyi2|] . It has also been suggested that 
this form is appropriate for localized electrons in semi- 
conductors 1^. In the present context ( p2|) is of inter- 
est because if the rotation matrix IR(i) were indepen- 
dent of time, our local rotation scheme would eliminate 
anisotropy to all orders, rather than just to first order 
in spin-orbit coupling. It is therefore natural to ask to 
what degree the fact that IR(t) depends on time spoils 
this hidden symmetry. 



For the particular form of anisotropic exchange in (22), 
the symmetric anisotropy term is, to second order in /3(i), 

TTabit) = - {PitfSab ~ Pa{t)l3b{t)) /2 + 0(/3(t)^). (23) 

For this ir(t), if we continue to take the pulse form ( p^ ) 
and (3(t) from ( p^ ) then the expression (|l^) can be eval- 
uated analytically. After performing the local rotation 



6Acot(A/2) - 12) 
X {Pl5ab-Piafiib)+0{(it). (24) 



—I _ 8Jo(A)2 



Thus even for the rotated exchange (|22|), if the rotation 
depends on time we are still left with residual second- 
order anisotropy after a pulse. 

To summarize, we have studied the effects of 
anisotropic corrections due to spin-orbit coupling on 
quantum gates produced by pulsing the exchange interac- 
tion between two spin-1/2 qubits. These quantum gates 
are parametrized by an effective Hamiltonian that com- 
mutes with itself at different times and produces the same 
quantum gate as a given pulse. Expressions for the vari- 
ous parameters in this effective Hamiltonian are obtained 
perturbatively in powers of spin-orbit coupling and used 
to shape pulses that effectively eliminate first-order spin- 
orbit corrections to quantum gates. The ability to reduce 
spin-orbit effects from first order to second order should 
be useful for any quantum computing scheme which relies 
on isotropic exchange. 
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